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Abstract 

We compute numerically the distribution of energies Q(E, N) for the XY-model 
with short-range and long-range interactions. We find that in both cases the dis- 
tribution can be fitted to the functional form: Q(E,N) ~ exp(N<j)(E,N)), with 
(f>(E, N) an intensive function of the energy. 
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In 1988 C. Tsallis introduced a generalized g-entropy given by [1] 

V w n q - 1 
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where pi, (i = 1, . . . , W) is the probability of the microscopic configuration i, 
and the parameter q characterizes the degree of non-extensivity of S q !^\ Non- 
extensive systems with long-range interparticle interactions are good candi- 
dates to be studied under the generalized thermostatistics derived from this 
g-entropy [2]. One such system which has received special attention recently 
[3-5] is the XY model with long-range interparticle interactions. This model 
is described by the Hamiltonian 



1 A quantity is said to be non-extensive if its asymptotic dependence on the sys- 
tem's size iV is non-linear. On the contrary, extensive quantities exhibit a linear 
dependence on N. 
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where {9i}i=i t ...,N are angle-type variables, and the sum runs over all distinct 
pairs of sites on a d— dimensional regular lattice of lineal size L = N 1 ^ with 
periodic boundary conditions. The distance between the sites i and j is r^, 
and the parameter a sets the interaction range. The short-range interaction 
XY model is recovered in the limit a = oo. To determine the range of param- 
eters leading to non-extensive features, it is instructive to consider the scaling 
behaviour of the mean energy per particle [6], 

Z~N = l + d'jdrr^r-"= Nl - a ' d - a J i (3) 
N J l-a/d K ' 

Therefore, for a > d (including the limiting case a = oo) the energy scales as 
E ~ N and the system behaves linearly, whereas in the non-extensive regime 
a < d, the energy scales as E ~ JSf 2 ^ a / d , i n ^he limiting case a = d the 
energy per site scales as the logarithm of the system size: E ~ N\nN. The 
number of microscopic configurations with energy between E and E + 5E is 
equal to Q(E, N)5E, where Q(E, N) is the density of states. Our aim here is 
to compute numerically Q(E, N) in order to study its asymptotic dependence 
on N. 



Recent numerical studies of a dynamical version of the long-range XY model 
(endowed with appropriate kinetic energy terms) have shown the existence of 
metaestable states exhibiting a Tsallis' maximum entropy distribution of ve- 
locities [5]. The lifetime of these metastable states increases with the system 
size, thus suggesting that it diverges in the thermodynamic limit. Unfortu- 
nately the relaxation time needed to arrive to these metastable states also 
increases with the system size, making it very difficult to study numerically 
the proposed non-commutativity: linijv^oo lim^oo 7^ lim^oo lim^v^oo which is 
expected to appear in the thermodynamic limit in these kind of non-extensive 
systems. 

Here we are going to consider an angle of this subject involving a recently 
conjectured relation [7] between (i) the range of applicability of Tsallis Statis- 
tics and (ii) the functional form of the density of states Q(E,N). For most 
systems studied so far one has Q(E,N) ~ exp(N(j)(E, N)), with <fi(E,N) an 
intensive function, implying that the natural expression for an extensive en- 
tropy is the well known Boltzmann-Gibbs entropic form S = ln(Q(E,N)). 
However, things might be different if, instead, the behaviour of Q(E, N) were 
as Q(E, N) ~ (p{E, iV)iV M (again with ip(E, N) an intensive function). In such 
a case it is conjectured that the relevant entropic form is the one determined 
by Eq.(l), S q = (Q(E, Nf~ q -I) /{l-q), with fi(l-q) = 1. Our intention is to 
present briefly some results about the form of Q(E, N) for the XY model with 
long-range interactions. It is shown that an unrestricted counting of states 
leads to a functional form Q(E,N) ~ exp(N<f)(E)). Therefore, if ergodicity 
is satisfied, the Boltzmann-Gibbs entropic form is the appropriate one. This 
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suggests, in accord with the results reported in [5], that transient regimes 
(during which the system is not able to explore the whole range of energy 
values) may constitute an appropriate field of application of the nonextensive 
thermostatistics (see also [2]). The number of states Q(E,N) has been nu- 
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Fig. 1. Scaling relations of the number of states Q(E,N) of the one dimensional 
XY model, Eq.(2), for different values of the parameter a and N = 50, 100, 200, 400. 

merically computed using the so-called Histogram by Overlapping Windows 
(HOW) method [8-10]. Briefly, the HOW method works by confining the en- 
ergy histogram generation to a suitably small energy interval [Ei, E,i + AE], 
such that the values of the number of states for each energy in that inter- 
val are comparable. The process is repeated varying Ei until the full energy 
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range is covered. We have also applied a method recently developed by Wang 
and Landau [11] which samples the function Q(E,N) autoconsistently. Both 
methods give the same results within errors. 

We have obtained numerically Q(E, N) for the one-dimensional XY model 
defined in Eq.(2), for different system sizes N = 50, 100,200,400. Our results 
are shown in Fig.(l) in three cases: (i) the extensive, short-range limit a = oo, 
(ii) the non-extensive case with a = 0.8, and (iii) the non-extensive case with 
infinite range limit a = 0. It is clear from Fig.(l) that in all the cases con- 
sidered the generic behavior is Q(E,N) ~ exp(N<ft(E, N)) with an intensive 
function <j)(E,N) = <f)(E/NN) in which the energy appears rescaled accord- 
ing to Eq.(3). In this case, and according to the above conjecture, it appears 
that the equilibrium properties of this system would be described by the stan- 
dard canonical ensemble while the results of [5] apply to the metastable states 
developed during the transient dynamics. 

Finally, we note that the scaling functions (f)(x) appears to be independent of 
a in the range < a < d. This equivalence for this range of values of a was 
also observed with other properties of the model in [3,4] where it was shown 
that a mean-field description holds in this case. 
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